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Abstract 

We study the effects of thermal fluctuations on elastic rings. Analytical expressions 
are derived for correlation functions of Euler angles, mean square distance between 
points on the ring contour, radius of gyration, and probability distribution of writhe 
fluctuations. Since fluctuation amplitudes diverge in the limit of vanishing twist rigid- 
ity, twist elasticity is essential for the description of fluctuating rings. We discover a 
crossover from a small scale regime in which the filament behaves as a straight rod, 
to a large scale regime in which spontaneous curvature is important and twist rigid- 
ity affects the spatial configurations of the ring. The fluctuation-dissipation relation 
between correlation functions of Euler angles and response functions, is used to study 
the deformation of the ring by external forces. The effects of inertia and dissipation 
on the relaxation of temporal correlations of writhe fluctuations, are analyzed using 
Langevin dynamics. 
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1. Introduction 

Small circular loops of DNA (plasmids) play an important role in biological processes such 
as gene transfer between bacteria and in biotechnological applications where they are used 
as vectors for DNA cloning . The simplest minimal model that captures both the topology 
and the physical properties of such an object is that of an elastic ring that has both bending 
and twist moduli. This model was used in a recent study of writhe instability of a twisted 
ring [0, ^ . However, since this work focused on the mechanical aspects of the problem 
and did not consider the effects of thermal fluctuations, it cannot be directly applied to 
plasmids and other microscopic rings. The consideration of fluctuations is important since 
they dominate the physics of macromolecules and determine their statistical properties, 
such as characteristic dimensions, djTiamics in solution [Q, kinetics of loop formation and 
dissociation of short DNA segments and molecular beacons 0. Recently, we developed 
a theory of fluctuating elastic filaments, with arbitrary spontaneous curvature, torsion and 
twist in their stress free state [0. Since topological constraints were not taken into account 
in this work, our analysis was limited to open filaments and could not be directly applied 
to the study of closed objects that have the topology of a ring. 

The present paper is an expanded version of a letter in which we presented the solution of 
this problem for weakly fluctuating rings |^ . The analysis of reference |^ is generalized to the 
case of ribbonlike filaments, with two principal axes of inertia in the cross sectional plane. We 
calculate the correlation functions of Euler angles, and use them to obtain other statistical 
properties of fluctuating rings, such as mean square spatial distance between points on the 
ring contour, and radius of gyration. Analytical expressions for the complete probability 
distribution function of writhe fluctuations and for all its moments, are derived. A crossover 
length scale is found, below which straight rod behavior dominates and the twist of the cross 
section with respect to the centerline is uncorrelated with the conformation of the centerline. 
Above this length scale the nonvanishing spontaneous curvature of the ring begins to play 
a role and twist rigidity affects the three-dimensional conformation of the centerline of the 
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ring. The correlation functions of Euler angles are used to predict the mechanical response 
to external torques and forces, and to examine the effect of spontaneous orientation of the 
cross section, on the deformation of ribbonlike rings. The dynamic correlation function of 
writhe fluctuations is calculated in both the inertial and the dissipative regimes. In the 
former case oscillatory decay of the correlations with time is observed. When inertia is 
negligible, the relaxation is monotonic and there is a transition from a short time regime in 
which the relaxation rate depends only on the bending rigidity, to a long time regime where 
the decay is affected by both bending and twist modes. 

In section ^ we present the generalized Frenet equations that describe the conformation of 
a filament, and introduce the elastic energy that governs its fluctuations about the stress free 
state. We express the curvature and torsion parameters that characterize this conformation, 
in terms of the Euler angles, and write down the elastic energy as a quadratic form in 
the deviations of these angles from their values in the undeformed ring. The topological 
constraints corresponding to a ring are introduced as integral conditions on the fluctuations 
of the Euler angles, and result in vanishing contribution of some of the lowest Fourier modes 
to the fluctuation spectrum. In section |^ we diagonalize the elastic energy, obtain the 
spectrum of normal modes and discuss their physical meaning. In section ^ we use this 
eigenmode expansion to calculate the correlation functions of Euler angles. We study the 
dependence of the correlators on physical parameters such as bending and twist rigidities, 
and on the spontaneous orientation of the principal axes of inertia of the cross section with 
respect to the plane of the ring, and discuss the geometry of typical configurations of the 
ring. In section |^ we derive explicit expressions for the orientational correlation function of 
the tangents to the ring, root mean square (rms) distance between points on the ring contour, 
and radius of gyration, in terms of correlation functions of Euler angles computed in the 
preceding sections. In section ^ we express the writhe and twist numbers that characterize 
an instantaneous configuration of the ring, in terms of Euler angles. We then use the 
correlation functions of Euler angles to calculate the probability distribution function of 
writhe fluctuations and study its dependence on the bending and twist rigidities. We find 
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that the amphtude of writhe fluctuations exhibits a crossover from small scale, straight-rod- 
like regime in which twist of the cross section has no effect on the spatial conformations of 
the centerline, to a large scale regime in which the two types of fluctuations become strongly 
coupled due to the spontaneous curvature of the ring. In section |^ we use the fluctuation- 
dissipation theorem that relates the previously calculated equilibrium correlation functions 
of Euler angles to the response functions, in order to study the linear response of a ring to 
small externally applied forces and moments. We show that the deformation of a ribbonlike 
ring depends in an essential way on the orientation of its cross section in the undeformed 
reference state. In section |] we derive the Langevin equations that describe both the inertial 
and the dissipative dynamics of Euler angles, and use them to study the effects of bending 
and twist rigidities and of the orientation of the cross section of the ribbon, on the frequency 
spectrum and temporal relaxation of its writhe modes. Details of the derivation of the 
Langevin equations and the calculation of the dynamic correlation functions, are given in 
appendices ^ and ^, respectively. In section |] we summarize our main results and discuss 
the domain of validity of our theory. 



2. General approach 

The general theory of fluctuating noninteracting elastic fllaments was presented in Ref. [0. 
To each point s one attaches a triad of unit vectors {tj(s)} where t^ls) is the tangent vector 
to the curve at s, and the vectors ti(s) and t2(s) are directed along the axes of symmetry 
of the (in general, noncircular) cross section. The spatial conformation x(s) of the fllament 
is given by the generalized Frenet equations 

= - J2jk ^ijkUJjtk, (1) 

together with the inextensibility condition, 

diL/ds = tg, (2) 
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where Cijk is the antisymmetric unit tensor and the parameters {ujj{s)} characterize the cur- 
vature, torsion and twist of the filament. The components of these vectors can be expressed 



in terms of the Euler angles 6, (p and ip: 

I 



ti 



cos Q cos (/? cos — sin sin 
cos Q sin cos ^/^ + cos sin ^/^ 
— sin Q cos 



(3) 



^ — cos Q cos ip sin '\\) — sin p cos ^ 



cos d sin sin t/; + cos cos 
sin ^ sin ■?/' 



sm cos p 



(4) 



V 



(5) 



Substituting Eqs. 



sin Q sin p 

cos 

@ into Eq. (|l]), the Frenet equations can be rewritten in the form: 



— = cui sm + UJ2 cosip, 
as 



— — smfc' 
as 



-oji cos '?/' + ci;2 sin '?/', 



(6) 



dilj 
ds 



sin 9 = {uji cos ip — uj2 sin ?/') cos ^ + tus sin 6. 



Solving these equations with respect to {tUj} yields 

dip , ^ , dO . ^ 
uJi = — — sm & cos ip + — smip, 
as ds 

df . n ■ > dO 
002 = —r o smip + — cos ip^ 
ds ds 

dijj ^dp> 
= — + cost^ — . 
ds ds 



(7) 
(8) 
(9) 



We assume that the centerline of the undeformed ring forms a circle of radius r in the 
xy plane, and that its cross section is rotated by angle iI)q{s) around this centerline. The 
Euler angles that describe this configuration are 



= 7r/2, = s/r, ipQ = ks/2r + ipoo, 



(10) 



where k is an integer and ipoo is a constant, independent of s. Eqs. (7) - (9) can be rewritten 
in the form 



t^oi = - (1/?") cosV^o, (^02 = {l/r)smtjjo and cuos = dtjjo/ds. 



(11) 



Although, in general, the stress free state of the ring can be arbitrarily twisted (e.g., because 
of intrinsic tendency of the filament to twist), in this work we will not consider spontaneous 
twist (cuos = 0), and taking k = we set ipo = i'oo (for brevity, we will denote this constant 
by ipo in the following). This angle characterizes the orientation of the principal axes of the 
cross section with respect to the plane of the undeformed ring. In the case of a circular cross 
section, all physical observables are independent of ipo and it is convenient to set ipo = 0- 
The corresponding Euler parametrization of the triad vectors is 



^ — sin(s/r) simjjQ ^ 



^01 



cos(s/r) sinipo 
— cos ipo 



^ — sin(s/r) cos tpo ^ 



^02 



V 



cos(s/rj 



^03 



J 



sin(s/r) 




(12) 



cos(s/r) cosipo 
sin ipQ 

In the absence of excluded-volume and other nonelastic interactions, the energy of a fila- 
ment is of purely elastic origin and can be represented as a quadratic form in the deviations 

6uJk = ujk - ujQk 

keT r^^r 



(13) 



where fc^ is the Boltzmann constant, T is the temperature, and the bare persistence lengths 
Ofc represent the rigidity with respect to the corresponding deformation modes. The above 
expression for the energy is based on the linear theory of elasticity and applies to defor- 
mations whose characteristic length scale (e.g., radius of curvature) is much larger than 
the diameter of the filament Since the persistence lengths are determined by material 
properties on length scales of the order of this diameter, they are the same as those of a 
straight rod. We conclude that ai and 02 are associated with the bending rigidities of the 
filament with respect to the two principal axes of inertia Ji and I2 (they differ if the cross 
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section is not circular), and that 03 is associated with twist rigidity. In the special case of 
incompressible isotropic rods with shear modulus fi, the theory of elasticity yields 

ai = Sfili/ksT, 02 = 3/i/2 //cfiT, and az = C/kBT (14) 

where the torsional rigidity C is also proportional to /i and depends on the geometry of the 
cross section (for an elliptical cross section with semi-axes di and d2-, C = Tx^d\d\/ {d\ + d'^). 
In this paper we will treat Oj as given material parameters of the ring. 

In the following we consider only small fluctuations of the Euler angles about their values 
in the undeformed state, Eq. (|l3). This approximation remains valid as long as the bare 
persistence lengths are much larger than the radius of the ring, i.e., ^ r. Expanding Eqs. 
(7) - (9) in small deviations from the stress free state, we find 

(5^ dde\ . , dSif 

duJi= { h — — smijQ — cost/^o, 

\ r ds J ds 

^ ( 6ip d69\ , d6(p . , 

'^^2 = I — + I cos?/^o + -^sm^o, (15) 

dStp 66 

duj-i = -3 • 

ds r 

It is instructive to relate the above parameters to the curvature k and torsion r familiar 



from differential geometry of space curves [jTOl- A circular planar ring has Kq = 1/r and 



To = 0. Expanding in small deviations about these values yields 



d5(p , , 66 d^du 

OK = — — and or = r = r (16) 

ds r ds^ 



As expected, fluctuations of the curvature represent bending deformations in the plane of the 
ring, and depend only on the angle (p that describes the rotation of the tangent to the ring, 
in the xy plane (see Eq. (^). Torsion describes deviations of the fllament from this plane, 
and its fluctuations depend only on the deviations of the angle 6 from tt/2. The speciflcation 
of the local curvature and torsion completely determines the conflguration of the centerline 
of any curved filament, and the Euler angle ip complements the description by specifying the 
rotation of the cross section about this centerline. However, the elastic energy can not be 



factorized into a sum of contributions due to deformation of the centerline and rotation about 
it. As will be shown in Section ^, ujs{s) defines the rate of twist and therefore the persistence 
length as is associated with twist. Twist represents not only rotation about the centerline 
(the dtp/ds term in Eq. (9)), but also contains a contribution due to the curvature of the 
centerline (the cos 6d(f/ds term in the above equation). Similarly, although inspection of Eq. 
(H) suggests that uji{s) and u!2{s) completely determine the variation of the tangent t3(s) as 
one moves along the contour, this variation depends on the main axes of the cross section at s 
(the vectors ti(s) and t2(s)), that themselves rotate with the cross section. This explains the 
■j/'-dependence of ui and UJ2 in Eqs. (7) and (8). The relation between the two descriptions 
{{6,(p,ip} and {uJi}) is a special case of the more general relation between Eulerian and 
Lagrangian descriptions in the theory of elasticity While the Euler angles describe the 



orientation of the triad {tj(s)} in the laboratory frame, the parameters uji{s) describe the 
local variation of this orientation as one moves along the curve, in the frame associated with 
the triad itself. The simple form of the energy, Eq. (|13[), is a direct consequence of this 
Lagrangian description. 

Substituting Eqs. (15) into the elastic energy, Eq. (|TBp, yields 



r2'n-r 

U = kBT ds 
Jo 



Y\d7^V) ^Y\d7) 

2' 



^ d69 d6^ as d6ip 69 

^ \ ds r J ds 2 \ ds r 

where the coefficients Ai are defined as, 



(17) 



Ai = ai cos^ ^0 + a2 sin^ -ipo, A2 = ai sin^ -ipo + a2 cos^ tpo, 

A3 = (a2 - ai) cos'?/'osin?/)o. (18) 

For a2 > ai, the constant Euler angle tpQ measures the angle between the major axis of 
inertia and the xy plane. The case 1^0 = ^ (V'o = ^/2) corresponds to major axis that lies in 
the xy plane (normal to the xy plane). The coefficients Ai obey the relations 

A1A2 -Al = aiaa, Ai + A2 = + ag (19) 



The periodic boundary conditions on the Euler angles 



S9{2TTr) = 5e{0), Sijj (27rr) = (5V^(0), 5ip{2TTr) = 5(p{0) 



(20) 



are supplemented by the condition that the ring is closed in three dimensional space, 
x(27rr) = x(0). Using Eq. this condition can be recast into an integral form, 



/ dsdt^is) = 0. 
Jo 



(21) 



For small deviations from equilibrium we get from Eq. (^), 



^ —5(p{s) sin(s/r) ^ 



6ts{s) 



6(p{s) cos(s/r) 
-6e{s) 

and the boundary conditions can be written as 



V 



(22) 



2nr 



dsSe{s) 



2irr p2TTr 

ds5ip{s) cos{s/r) = / ds6(p{s) sm{s/r) = 0. 
Jo 



(23) 



Since the deviations of the Euler angles are periodic functions of s, they can be expanded 
in Fourier series 



(24) 



for each rj = 6,ip, ip, where the sum goes over all positive and negative integers n. The 
boundary conditions, Eqs. (^), can be expressed as conditions on the Fourier coefficients. 



m = (^(1) = 0. 



(25) 



Substituting Eqs. (|2|) into Eq. (17) we find. 



U 



27irkBT 



^(0) +(Ai + a3) 19(1) +4,(1) 



+ 



(26) 



n=2 



^ IAi in6{n) + tlj{n) + A2n'^ \ip{n)\ 



2A^ ^in6{n) + ilj{n) inip{—n) + 03 inil){n) — 6{n) 



The energy does not depend on modes ■0(1) = —^^(1) and <^(0) that correspond to rigid 
body rotation of the entire ring, with respect to axes lying in the plane of the ring and 
normal to it, respectively. 

The quadratic form inside the sum over n in Eq. (26) can be represented as a matrix in 
the space spanned by the Fourier components 9{n), (p{n) and %l]{n) (this applies to n > 1; 
the cases n = 0, ±1 will be considered separately). 



Q(n) = 



^ Aiin? + as A^in? —i {Ai + as) n ^ 

A^-n? A2n^ —iAsn 
^ i {Ai + as) n lA^n asu^ + Ai ^ 



(27) 



3. Spectrum of fluctuations 



In order to obtain the spectrum of fluctuations of the ring, we diagonalize the free energy, 
Eq. (26), by expanding the Fourier components fj{n) in the eigenvectors r]k{n) of the matrix 
Q(n), 



(28) 



where rj — 0,(f,ip and r)k{n) is the 77-th component of the eigenvector rjkin) — 
{9k{n), ipki'n),ipkin)} of the quadratic form, Eq. (26), corresponding to the eigenvalue Xk{n). 
They are normalized by the conditions. 



Y.„^k{n)r]i (-n) = 6ki, I].%(r^)% (-n) = 5^^/. 



(29) 



Expanding the elastic energy in the eigenmodes gives 



U = 



^k{n) \ck{n)\ . 



(30) 



' n=0 k 

The three eigenvalues Afc(n) corresponding to the Fourier mode n, are the roots of the 
characteristic cubic polynomial. 



- 62A^ + 61A - 60 = 0, 



(31) 
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with coefficients 



bi = {aia2 + + ^1^3) (n"^ ~^ ^) ~ ^^^^ ("^""^ ~ ^) ' 
&2 = (cti + 0-2 + as) n"^ + Ai + 03, 

where we used Eqs. (18) and (19) to simpUfy cumbersome mathematical expressions. Since 



the matrix Q(n) is Hermitian, its eigenvalues are real. 



Inspection of Eqs. (|3l|) and (32) shows that Xk{—n) = Xk{n) and that all eigenvalues 



with n > 1 are positive. Because of the boundary conditions, Eqs. (25), there are only two 
independent normal modes corresponding to each of the cases, n = and n = 1. In order 
to understand the physical meaning of these modes, we introduce the components of the 
Fourier transforms of the curvature and torsion, Eq. (|16|), 



in , , _ , , n'^ — 1 ~ 



Kin] 



r 



(p{n) and r(r?,) = 6{n). (33) 



Substituting the boundary conditions 9{0) = = into the above expressions we con- 
clude that for modes with n = and 1, both 6k, and 6t vanish and, therefore, these modes 
do not affect the planar circular configuration of the centerline of the ring. There are two 
zero energy modes that correspond to symmetry operations on the undeformed ring. One 
n = mode, with eigenfunction v^i(O) = 1; ipi{0) = 0, describes the rotation of the ring 
about the z axis. One n = 1 mode, with eigenfunction 6*1(1) = 1; ipi{l) = —i, corresponds 
to rotation of the ring about an axis in the xy plane. The two remaining modes have an 
energy gap and are twist modes that leave the centerline undisturbed. The n = mode with 
eigenfunction (p2{0) = 0; 'ip2{0) = 1 has an eigenvalue A2(0) = Ai, and describes uniform 
twist of the ring about its centerline. Since this eigenvalue does not vanish for arbitrary 
ai and 02, we conclude that uniform twist of a ring costs energy even if the ring has a 



circular cross section. This conclusion agrees with reference []T2[ , where the dynamics of the 
uniform twist mode was studied. The n = 1 mode with eigenfunction ^2(1) = 1; V^2(l) = i 
has the eigenvalue A2(l) = 2(y4i + a^) and corresponds to rotation of the ring with respect 
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to an axis that passes through the centerhne and hes in the xy plane, accompanied by twist 
of the cross section by the angle ip that varies periodically (as cos(s/r)) along the contour 
of the ring. The dynamics of this mode was studied in reference |]T3| . 

In the limit n ^ 1, fluctuations of the three Euler angles are decoupled and \k{n) ^ akU^. 
In general, each normal mode of the ring corresponds to fluctuations of all three Euler angles, 
5ip{s) and 5il){s), and describes a complex three-dimensional configuration. 

The eigenvalue problem is simplified for a circular cross section (02 = cti), or when 
the cross section is asymmetric but ipQ = Q (the case z/'o = 7r/2 is reduced to ipo = Q by 
the substitution ai ^ 02). In these cases the mode 5(p{s) decouples from the other two 
modes and has the spectrum \i{n) = aiTi^ {n 7^ ±1)- This mode corresponds to bending 
fluctuations that lie entirely in the plane of the ring. The other two modes are linear 
combinations of 66{s) and Sip{s), with eigenvalues 



A.,3(n) = ^^(n^ + 1) ± y (^^) (n^ + 1)^ + ^n^a.a,. (34) 

Eq. (^4|) can be further simplified in the limit of large rigidity with respect to twist, 03 ^ 02, 
in which case 



n^ + l '~ ^' (35) 

Xzin) = as{n'^ + 1) for n >2. 



\2{n) = a2 2 I 1 for n > 1 



In the opposite limit as <^ 02, the eigenvalues can be found by substituting 02 as in Eq. 

Inspection of Eq. ( [5^ ) shows that Xsin) vanishes identically for all n when 03 = (this 
statement applies even to rings with noncircular cross section - see Eqs. (^Tj) and (32)), 
indicating that the amplitudes of the corresponding fluctuation modes grow without limit 
in the absence of twist rigidity. Examining the expression for the elastic energy, Eq. (17), 
we conclude that these zero energy modes correspond to fluctuations for which 

d6e/ds = (36) 

In the absence of twist rigidity, twist fluctuations carry no energy penalty and the angle of 
twist of the cross section {dif)) can always adjust itself to arbitrary deviation of the center line 
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from the plane of the unperturbed ring (69), so that this condition, Eq. (|36|), its satisfied. 
The presence of an infinite number of zero energy modes means that twist rigidity {a-^ ^ 0) 
is absolutely essential for stabilizing the ring against out-of-plane fluctuations, and that 
bending elasticity alone can not suppress this instability. 



4. Correlations of Euler angles 



Applying the equipartition theorem to Eq. (pOD, we get 



(cfc(n)cfe'(-n')) 



(37) 



Using expansion p8| ) and averaging with the help of Eq. (|37D , the correlation functions of 
Euler angles can be expressed in terms of the eigenvalues Afc(n) and the eigenfunctions rik{n) 
of the Q(n) matrix: 



{6r] {s) 6r]' {s')) = ^ e^"(-^')/r (^(^)^'(_^)) = L ^ ^Ms^s>yr ^ Vkin)r]',{ n) 



(38) 



where 67], 61]' = 66, 6ip, dip. Care should be exercised in evaluating the above expression, 
when considering the contribution of the modes with n = 0, ±1, since modes with vanishing 
eigenvalues should be excluded. A straightforward calculation gives 

^ ^in{s-s')/rY^ Vk{n)rik{-n) 



n=0,±l 
1 



^0 0^ 




1 



Afc(n) 



1 



/ (s-s'\ . (s-s'\ \ 
cos — sm 



^1 + 03 



r J 








V 



s — S \ I s — s 

sm cos 



r 



(39) 



r y y 

where rjkfjk denotes the direct product of two vectors rik, the t]-)]' component of which is rjkri'i^. 
For 7^ 0, ±1 we flnd 



(40) 
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where Q ^(n) is the inverse of the matrix Q defined in Eq. 



Q-\n) 



^ I If —^3 

a3{n^ — 1)^ a±{n'^ — l)'^ aia2(n^ — 1) 



(- + - 



a±J (n2-l)2 



V as a±J 



iA-> 



aia2(n^ — 1) 



(41) 



3 nz I i 

(n^ — 1)2 aia2n(n2_i) a3(n2 — 1)2 a^(n2 — 1)2 y 



and where 



1 



cos^ ipo sin^ ipQ 
a± ai 02 



1 sin^ ipQ ^ cos^ ipo 



(42) 



a|| ai 02 

Effective persistence lengths 03 and a± control both fluctuations perpendicular to the plane 
of the ring and fluctuations of the twist angle ip, and an controls fluctuations in the plane 



of the ring. Using Eqs. (|38|) and (PD, we obtain the correlation functions of Euler angles 
(here s = |s2 — Si| , < s < 27rr) 



{69{s,)69 {S2)) 

(5^(S1)5^(S2)) 



r cos(s/r) 

— . + ■ . 

n Ai + a-i TTtts 



■/a - 



A - +— /i - 



7ra_L 



Trail 



+ - 



r cos (s/r) r 



+ 



r sin (2'?/'o) / I 1 



27T 



a2 tti 



/i - + — h - 



7ra_L 



(43) 



r sin(s/r) ^ r ^ 1 ^ 1 



TT Ai + as TT Vas a_L 
r sin (2?/'o) 



/i(r)-/3(- 



2n 



/5 - 

a2 ai 



where we deflned, for < a; < 27r, 



f2{x) 

fsix) 



^ n cos nx 
h in' - if 
cosnx 



(tt — xY 



24 ^ 16 



TT — X 



cosx 



sm X, 



E 

n=2 

CO 



/ n2 2 

TT — X TT 



cosnx 



(vr — x) 



12 

2 2 



cosx, 



TT- _ 3 

24 ~ 16 



n sm nx 
(n2 - 1)2 
sin nx 



E 

?i=2 

00 

S^(n2-1) 



(vr — x)^ TT^ 



1 



smx 



24 16 

vr — X 



vr — X 1 

cosx H smx , 

4 2' 



smx, 



(44) 



cosx — 1) . 
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Inspection of Eqs. (43) shows that the bare persistence length associated with the twist 
rigidity, 03, plays a fundamentally important role: fluctuations of the angles ip and 6 and the 
correlation between these angles, diverge in the limit 03 0! Therefore, simplified models 
of elastic filaments with nonvanishing spontaneous curvature that do not take into account 
twist rigidity, can not describe fluctuations and elastic response of the ring. This is not 
the case for a straight rod, whose spatial fluctuations can be successfully described by the 
wormlike chain model Jill ("with 03 = 0). The reason for the difference stems from the fact 
that the elastic energy of straight rods contains no coupling between the angles that describe 
the spatial conformation of the centerline {9 and ^p) and the angle that describes the twist 
of the cross section about this centerline {ip). When twist rigidity vanishes (03 = 0) there 
is no energy penalty for twisting the cross section about the centerline and the amplitude 
of twist fluctuations of the cross section about the centerline diverges, but the presence of 
bending rigidity (ai, 02 7^ 0) suffices to suppress spatial fluctuations of the centerline about 
its straight stress free configuration. For rings, the elastic energy in Eq. (17) contains cross 
terms in the angles dip and 59 that couple both types of fluctuations. Inspection of Eq. (17) 
shows that when 03 = 0, fluctuations with d59/ds + 5ip/r = have zero energy cost (see 
Eq. (PBD) and, since in the absence of twist rigidity the angle Sip can always adjust itself 
to satisfy the condition 6ip = —r d69/ds, for 03 = there is no elastic energy penalty for 
out-of-plane fluctuations of the ring and the amplitude of such fluctuations diverges. We 
conclude that standard wormlike chain theories in which only bending rigidity is taken into 
account, can not model fluctuating rings. 

In Figs. 1 - 2 we plot correlation functions of Euler angles, for a ring with circularly 
symmetric cross section. Substituting ai = a2 in the expressions for the angular correlators 
in Eqs. (43) we find {69{si)6{p {32}) = {6Lp{si)6ip {32}) = 0. The physical reason for this 
behavior becomes clear when one recalls the discussion of the eigenvalue problem for a 
ring with circularly symmetric cross section (see Eq. (0)). In this case, fiuctuations of 
(p{s) decouple from those of the other two angles and therefore, cross correlation functions 
involving 6(f vanish identically. In Fig. 1 we consider the case oi = 02 -C 03, i.e., twist 
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rigidity is much larger than the that of the bending modes. The diagonal angular correlation 
functions are oscillatory functions of the contour distance, with maxima at |s2 — Si| = 0, 
vrr and 27rr (they are symmetric with respect to reflection about the point |s2 — si| = vrr). 
These behaviors result from interference of two wave packets propagating along two opposite 
directions along the ring. As a consequence of the large twist rigidity, the correlator of the 
twist angle is always positive, while {69{si)69 (32)) and {6ip{si)6ip (32)) fluctuate around 
zero. The cross correlation function, {69{si)6ip{s2)) , vanishes as |s2 — si| — 0. The physical 
reason for this surprising behavior is that a short segment of the ring confined between these 
points can be considered as a nearly straight incompressible rod. Since twist of such a rod 
does not produce any deformation, local fluctuations of twist and of the other two modes 
are not correlated with each other. For larger contour separations, spontaneous curvature 
begins to play a role and fluctuations of 6 and ip become coupled. This is a manifestation 
of the crossover from small scale (twist and spatial conformation fluctuate independently) 
to large scale (coupled twist and centerline fluctuations) behavior, that will be discussed in 
greater detail in section ^. 

In Fig. 2 we present the case of small twist rigidity, Oi = 02 S> 03. The twist cor- 
relation function develops four nodes (i.e., points at which it vanishes) and, at the same 
time, its amplitude is strongly enhanced. In Fig. 2 we did not plot the correlation function 
{6ip{si)6ip (S2)), since it depends only on the bending rigidities (see the second of Eqs. (43)) 
and is therefore the same as in Fig. 1. Fig. 3 deals with the case of an asymmetric cross 
section (or asymmetric rigidity in the cross sectional plane), ai 7^ 02. The cross correlations 
{S9{si)6ip {S2)) and {Sip{si)S4' (32)) no longer vanish (for ipQ 7^ 0,7r/2), even though their 
amplitude is much smaller than that of {66{s)6il!{0)). Since the arguments presented in the 
preceding paragraph apply here as well, the two cross correlation functions involving Sip 
vanish as S2 Si. The cross correlation function {66{si)6ip {S2)) behaves in a way similar 
to that of the diagonal correlation functions and is symmetric about |s2 — si| = vrr. 

We would like to comment on the physical meaning of fluctuations of the angle ip{s). We 
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find from Eq. (43) 

/ „v 9t 1 r / « \ " 

(45) 



' a|| TT a| 



1 - 2 cos ( - 

r 



wliere tlie "parallel" persistence length a\\ is defined in Eq. (0), and where s\\ = s(l — s/27rr) 
is the effective contour length for parallel connection of two segments, one of length s = 
\s2 — Si\ and the second of length 27rr — s (analogously to parallel connection of resistors in 
an electrical circuit). The effective elastic modulus between points Si and S2 is proportional 
to 

— = - + , or S|| =s(l . (46) 

S|| s 27rr — s \ 2'Kr J 

The second term on the rhs of Eq. (|45|) arises due to subtraction of the contribution of the 

mode <^(1) because of the closure of the ring. Eq. (|45|) describes the Brownian fluctuations 

of phase (y9(s) on a circle, with effective "diffusion" coefficient a^^^. This means that the 

angle can jump discontinuously from point to point and therefore, the amplitude of its 

derivative dip/ds diverges. Since d^p/ds is the local curvature of the filament (see Eq. (|16|)), 

we conclude that <^[5/t(s)]^^ ^ oo. A similar calculation for the second derivative of the 

angle 6 shows that its amplitude diverges and therefore (^[5r(s)]^^ — oo as well. The above 

divergences are eliminated by a cutoff on length scales of the order of the thickness of the 

filament and, on length scales larger than this diameter, the contour of the ring remains a 

smooth and continuous curve in the process of thermal fluctuations. 



5. Spatial correlations and radius of gyration 

We proceed to calculate the correlation function <^[x(si) — x (52)]^) that measures the mean 
square spatial separation between points si and S2 on the contour of the filament. Integrating 
Eq. (^, yields x(si) — x(s2) = t'i{s')ds' and we can express this correlation function in 
terms of the correlator of tangents to the ring, at two arbitrary points on the contour, 
(13(5') • t^{s")). We show below that this orientational correlation function of the tangent 
vectors can be expressed in terms of correlation functions of Euler angles. Expanding the 
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vector ts to second order in deviations of Euler angles, St], from their unperturbed values 
gives, 

1 



ta = to3 + Sts = 66t'Q^+6ipt'Q2 + 



1 - - [66' + 6^' 



^03, 



(47) 



where the vectors t'^-^s) are defined by 



toi(s) 



— sin(s/r) 



5 "^02 (■5) 



V 



cos(s/r) 




cos(s/r) 



to3(s) 



sin(s/r) 




(48) 



When ipo = 0, these vectors coincide with the vectors of unperturbed triad, Eq. (|T^). Using 
Eq. (^7]) we find (in matrix notation) 



(49) 



MsMs^)) = (1 - {69') - {6ip')) to3(si)to3(s2) 

+ {69is^)6ifis2))t',,is^)tUs2) + {6cpis^)69is2)) tUsi)tUs2) 



where tojtoj denotes the direct product of two vectors toi and toj. The correlation functions 
of the Euler angles that appear in the above expressions are given in Eq. (43). As expected, 
the normalization condition for unit vectors, (t3(s) •t3(s)) = 1, is satisfied up to terms of 
second order in 61]. 
Using the equahty 



JO 



(50) 



valid for any even function /(x), we obtain 



x(si) - X {S2)f) = r ds' r ds" {t,{s') . t,{s")) 

' Jsi Jsi 



2r' 



1 — cos I - 

r 



r 

TT 



—9\\ (-) + —9± (-) + —93 (- 

ail \r/ a± \r/ 03 \r 



(51) 
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where s=|s2 — < s < 27rr. and where a_|_^ is defined in Eq. (^21) . The functions g\\, 

g± and are given by 



X 



9\\{x) = 2 (x-u) [/2(0) - /s (u)] cosudu 



[1 + cosx) (ttx ~ ^ j ^ 2 '-''^ '^°^-^)^ + 2 (vr — x) sinx + 2, 



X 



g:i{x) = 2 (x-u) [/i(0) cosn - /i (n)] du 







1 / X^ \ 1.x 1 

= -- I xvr — ^ + 1 I cosx + - (tt - x) sinx + -, (52) 

93{x) = 2 f (x-u) [/3(0) COSM - fs (u)] du 





If a;2 \ 3 , , . x2 3 

= -- I xvr ^ + 3 1 cosx + - (tt - x) smx - xvr + — + -. 

For small x <C 1 we have g\\{x) ^ g±{x) — 7rx^/12 and (73 (x) ~ x^/32 ^ 5'||(2;)- Com- 
bining these expressions into Eq. (51), we conclude that the lowest order corrections to 
the straight line result, <^[x(s) — x(0)]^^ / ~ '^^P^^'^ ^^^^ '^^ effective bending 
persistence length 2/ (^ai^ + a^^^, in agreement with the wormlike chain model. For general 
s this correlator depends on all the bare persistence lengths, ai, 02 and a^. 

In Fig. 4 we plot the mean square distance between two points on the ring contour, 
<^[x(s) — x(0)]^^, as a function of s, in the interval < s < 27rr. As expected, it in- 
creases parabolically with s (straight rod behavior for small s) and exhibits a maximum at 
s = vrr (the maximum is determined by the geometry of the undeformed ring). Fluctuations 
suppress this maximum in a way that depends on the various rigidity parameters. Thus, 
decreasing the twist rigidity, 03, has a much smaller effect on the amplitude of the maxi- 
mum, than decreasing the bending rigidities ai or 02. The origin of this effect is that twist 
rigidity does not affect the spatial conformations of a short segment of the ring that can 
be considered as a nearly straight incompressible rod. Therefore, twist fluctuations affect 
only the conformations of long segments, for which deviations from a straight rod become 
signiflcant (compare solid line and boxes in Fig. 4). 
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The radius of gyration is defined as 

1 



2nr 



ds 



X s - 



27rr 



(53) 



Averaging this expression over fluctuations, we can express (Ryj in terms of the two-point 



correlation function 



Rt 



^ j dSi j dS2 ([X(si) - X {S2)f) . 



(54) 



Using Eqs. (51) and (52) we find 



r: 











VStt 





3 r 
47r a± 



7 n\ r 
An QJ as 



(55) 



All the corrections to the unperturbed result (r^) are negative, and we conclude that fluctu- 
ations make the ring more compact. Since our weak fluctuation approximation is only valid 
in the range ^ r, these fluctuation corrections are rather small. Because of the small 
coefficient in front of the r/a^ term, the effect of twist fluctuations on the radius of gyration 
is relatively weak, but fluctuations diverge and the expression for the radius of gyration 
becomes unphysical in the limit of vanishing twist rigidity, 03 ^ 0. 



6. Writhe fluctuations 

The twist number Tw associated with a conflguration of the ring can be expressed through 
the Euler angles, 

Tw = ^ju,{s)ds = i^/"^ + cos^g) ds, (56) 

where we used the definition of the rate of twist, ^3(5), about the tangent vector, in terms of 
the Euler angles, Eq. (9). In order to understand the physical meaning of u^, consider the 
variation of the triad vector ti (or 12) as one moves an infinitesimal contour distance ds along 
the centerline of the curved filament. The projection of the vector ti(s + ds) on the cross 
section at s (the plane normal to the tangent t3(s)), rotates by an angle ujs{s)ds compared 
to its original direction, ti(s). Inspection of Eq. (|56D shows that this rotation consists of 
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two contributions. The first term corresponds to the contribution of a straight filament of 
length ds (the normal planes at points s and s + ds remain parallel to each other), whose 
cross section is twisted around the centerline, by an angle dip. The second term, cos 6dip/ds, 
arises due to the curvature of the centerline; since the cross sections at points s and s + ds 
are, in general, tilted with respect to each other, the projection of ti(s + ds) on the cross 
section at s will rotate by cosOdip. Notice that because of the interplay of the two effects, a 
curved filament can have zero twist even if dip/ds 7^ 0. This effect will be demonstrated in 
Section ^ (see Fig. 6). 

In addition to the twist of the filament that is closely associated with the rotation of 
the cross section about a curved centerline, and can be defined both locally (the twist 
"density" ^^3(5)) and globally (Tw), one can introduce an integral characteristic of the 
spatial configuration of the centerline that reflects its tortuosity, known as writhe number. 
In order to express the writhe number Wr of a given configuration of the ring in terms of 



Euler angles, one usually begins with the Fuller equation for the writhe of a closed curve [15 



d 



Wr = — i — — - — ds, (57) 



1 /• (tosxta) • — (tos+ts J 
2^ J 1 +to3-ts 

where x and • denote respectively vector and scalar products. In the above expression we 



made use of the fact that the writhe number of a planar circular ring vanishes |T6[. The 
above expression is valid as long as |tos ■ ts(s)| < 1 in the denominator, for all points s on 
the contour of the ring. This condition is satisfied in our work, since we only consider small 
fluctuations about a planar undeformed ring that lies in the xy plane. 

A more physically transparent definition of writhe is based on the existence of a topo- 
logical invariant of a ring, called the linking number |17| Lk. The total rotation of the cross 



section as one moves around the contour of the ring is given by 2'KLk where the linking 
number 

Lk = Tw + Wr, (58) 

does not depend on the conformation of the ring and is therefore a conserved quantity. In 
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the absence of spontaneous twist both the twist and the writhe numbers of a planar circular 
ring vanish, and Lk = Lk^g = 0. In general, since Lk is a constant for a given topology, 
6Lk = 6Wr + 6Tw = 0, and expanding the integrand in Eq. ( |56D in small deviations of the 
Euler angles from their spontaneous values in the unperturbed ring, the deviations of writhe 
and twist can be expressed as 

SWr = -6Tw = — / —^69ds = V in^{n)9{-n), (59) 
zvr JO as " 



where we used § d6ip = and § 50ds = (see Eq. (p3|) ). The last equality in Eq. (|59[) 
was derived using Eq. (12^) . Notice that the integrand in Eq. ( ^Of ) depends on the product 
of 6(p and 66, and we conclude that writhe deviations vanish both when fluctuations of the 
angles are confined to the plane of the ring {66 = 0), and when they are normal to this plane 
(6^ = 0). 

Since (^ip{n)6{—n)'^ is an even function of n, multiplying by n and summing over all 
positive and negative integer values of n yields 



{6Wr) = in [^{n)6{-n)) = 0. (60) 
The dispersion of the writhe number is given by 



6Wr^) = Wr^{n,. 

6{n)6{-n)) ((^(n)^(-n)) - (ip{n)6{-n ' ^ 



(61) 



where we excluded modes with n = and n = ±1, because of the boundary conditions 
^^(0) = <^(1) = 0. Using Eqs. (43) for the correlation functions of Euler angles, we obtain 
the mean square amplitude of writhe fluctuations at wavelength r/n, 

^ ' 7r2aia2a3 {n^ -if ^ ' 

Notice that the amplitude of writhe fluctuations diverges at as = and we conclude that 
twist rigidity plays an essential role in stabilizing the contour of the ring against writhe 
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fluctuations. The origin of tliis divergence is tlie same as tliat of the correlator {S9{s)S9 (0)) 
in Eq. (43) and has been discussed following Eq. (44). 

For large wavevectors, \n\ ^ 1, the mean square amplitude of the n-th mode of writhe 
fluctuations depends only on the bending persistent lengths, ai and 02. The physical reason 
is that on sufficiently small scales, the filament behaves as a straight incompressible rod 
whose properties do not depend on the twist persistence length 03 (see reference [|l4l)- In 
the limit Ai ^ 03^^, the writhe-writhe correlation function for a straight rod takes the 
form, 

Wr\q) = (63) 



where we defined the wavevector q = 2Txn/r. Eq. ( |63D is valid for straight rods when 
27r/g ^ Ob, where the persistence length a?, of the rod is defined by 

i = ifl + lV (64) 

2 \ai 02/ 

The crossover to a long wavelength regime in which writhe modes become affected by twist 



rigidity, takes place at a length scale = fyO'h/O'?, and, therefore, such a regime exists for a 

ring of radius r only if afe/a3 < 1. As a consistency check, notice that the straight rod case 

follows from the above expression for by substituting r = 00, and since diverges in this 

limit, Eq. ( |55| ) applies throughout the entire range of parameters. 

Substituting Eq. (|^) back in Eq. (|6lD yields 

-1 1 \ r2 /I 11 \ 



' ^ Stt"^/ aia2 Vo Stt"^/ 01103 



Notice that {5Wr^) ~ r^, in agreement with the scaling estimates in reference |T^. Indeed, 



since writhe is a quadratic form of 5ip and 59 (see Eq. (|59D), each of which has typical 



fluctuations of yr/a (a is a characteristic persistence length), the characteristic amplitude 
of writhe fluctuations is 5Wr r /a. 

The entire probability distribution of writhe can also be computed. Beginning with the 
formal definition of this distribution 



P [5Wr) = ( 6 



6Wr — ^ inip{n)9{—n) 

n^O,±l 
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(66) 



and using the exponential representation of the (5-function, yields 

P {5Wr) = / '^e'^^^'' (exp [k n(p{n)e{- 

kSWr J-J 



°° dk 
— ( 

oo 27r 

°° dk 



-n] 



2ti 



det Q(n) 



n^o ii [Q(n) + A;nY] 
where the matrix Y is defined as 



'^0 -1 0^ 



(67) 



TT 



1 




(68) 



Calculating the corresponding determinants gives, 



P {5Wr) 



°° dk 



AkSWr 



(69) 



^^'i^Y{n=2[^ + Wr\n)k^\ 
This integral can be calculated by expanding the integrand into partial fractions and we get 



P{5Wr) = ^7r(n) 



n=2 



An) = n 



1 - 



2Wr{n 
Wr^{k) 



exp 



\6Wr\ 
Wr{n) 



n -1 



Evaluating the product in Eq. (ffTD yields, 

Tcn'^ {n^ + a) (1 + a) -Ja 



(70) 
(71) 



Txin) 



a\n] 



2 irfi — 1) sinh [t^^Jol) 



(72) 



(n2-2)Ai -ag 



Ax + n^az 

The above expression for -nin) can be used to calculate all the even moments of writhe 
fiuctuations (odd moments vanish due to the radial symmetry of the undeformed ring). 



mr^) = k\Y,'K{n)Wr^(ji), 7r(n) = 1, fc = 2,4. 



(73) 

n=2 n=2 

Since moments with A; > 2 do not vanish, it is obvious that the writhe distribution is 
not Gaussian. Furthermore, inspection of Eq. ( |7DD shows that the distribution has an 
exponential tail at large 5Wr. For strongly writhing rings, {5Wr'^)^^'^ ^ |5iyr|, the n = 2 
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term dominates the sum in Eq. ([70|) and the free energy F = — ksT In P (SWr) is given by 
(up to logarithmic corrections), 

6Wr ^ I ^^^^ 2\V2 



for UWr') < \mr\ < 1. (74) 



The second inequahty, |5H^r| <^ 1, follows from the assumption that the deviations of Euler 
angles from their equilibrium values, are small. 

The writhe distribution function can be written in the form 

P i6Wr) = ^-^p ( ^-^5Wr, ^\ . (75) 
r \ r as J 

Plots of the dimensionless function p{x, Ai/a^) for ^1/03 = 0.1, 5 and 20 are shown in Fig. 
5. As intuitively expected, the probability of large writhe fluctuations (large TCy/oio^dWr /r) 
decreases with increasing twist rigidity (decreasing Ai/a^), but the effect saturates for 
Ai/a^, < 0.1. The shape of the curves bears close resemblance to the results of recent 
computer simulations [|19| . 



7. Elastic response of the ring 

According to the fluctuation-dissipation theorem, the correlation functions of the Euler 
angles determine the elastic response of the ring to external distributed torque M(s), apphed 
along its contour. In the following, we use this information in order to study the deformation 
of the ring by external torques and forces. Since we are not interested in rigid body rotation, 
we assume that the total torque on the ring vanishes, i.e., 

^dsM(s) = 0. (76) 

The deviations of the Euler angles from their unperturbed values are given by 

^V{s) = ^Y.r,Jds'{6r^{s)67^'{s'))M,,{s'), (77) 

where t], t]' = 9,(p,ip and M^/ are the corresponding components of the torque. In order to 
calculate the elastic response to external force, F(s), applied to the centerline of the ring, 
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we rewrite the work done by this force, W = § dsF (s) • (5x(s), in the form 

W = f dsm{s) • 5t3 (s) , where m(s) = / ds'F{s'). (78) 



Since we are not interested in translation of the ring as a whole, we assume that the total 
force acting on the ring vanishes, § dsF{s) = 0, which means that the function m(s) is 
continuous at s = 0, i.e., m (27rr) = m(0) = 0. Using Eq. (j2^) we can recast the expression 
for the work done by the force, Eq. (^), in the form 

W = f ds {[—mi (s) sin (s/r) + m2(s) cos(s/r)] 6(p{s) — (s) S9{s)} . (79) 



Inspection of this equation shows that in the presence of external force we have to modify 
the expressions for the moments 

My,{s) M^{s) — nil (s) sin(s/r) + 1712(3) cos(s/r), 
Me{s)^Me{s)-m3{s), (80) 

in Eq. (^). The condition that the total torque due to the external force vanishes, Eq. 
([76|), imposes additional conditions on the force F(s). Upon some algebra, these conditions 
can be written in the form, 

/•2irr r2T!r 

Ft{s)ds = sF3{s)ds = 0, (81) 

where Ft{s) = F{s) ■ to3(s) is the tangential component of the force F. Inspection of Eq. 
(80) shows that a small force with mi{s) = —Fr{s)rcos{s/r), 7/12(5) = —Fr{s)rsm{s/r) 



and 1713(3) = does not deform the ring. From Eq. (|78[) we find that Fr{s) is the radial 
component of the force F(s), while its tangential component is Ft{s) = rd [Fr{s)] /ds. This 
tensile force balances the contribution of variation of the radial force Fr{s) along the contour 
of the ring and prevents buckling until a critical value of the radial force is reached. 

Equations ( ^Tf ) and ( ^8]) together with Eqs. ([77|), (80) and inextensibility condition, Eq. 



(0), determine the conformation of the deformed ring, under the action of small external 
torque and force. As an illustration, consider the deformation of the ring under external 
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forces F applied to two opposite points s = ti/2 and s = 37r/2 on the ring contour, 

Fi (s) = F[5{s- 7r/2) - 5 {s - 37r/2)] . (82) 

Using Eqs. (^), (^) and (80), we obtain the following expressions for the resulting varia- 
tions of Euler angles. 



271 \a2kBT aiksT J \r 



2 

r , f s 



6^{s) = F —-h^ - , (83) 

7ia\\kBT \rj 

sin 2^0 ^ 1 , (s 

S^{^s) = -F — — - — — ]h^{- 
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where 

he{x) = = E ^^^^sm(2M = --xcosx, 

~ (-1)^= TT 

h^x) = E ^^4^2 _ 1) sin(2A;x) = x - - sin x, (84) 

h^ix) = — 2 > ——7; — cos[2kx) = 1 cosx — -xsinx. 

^1 (4A;2 - 1)2 ^ ^ 4 4 

The above series are calculated in the interval |x| < 7r/2. Using the periodicity condition, 
hi{x + vr) = hi{x), the functions hi{x) can be extended outside this interval. Inspection of 



Eq. (|l^) shows that the persistence lengths Oj are inversely proportional to temperature. 
Since temperature enters Eq. (|83| ) only in combinations aiksT, it cancels from the above 
expressions and can affect the results only through temperature dependence of elastic moduli 
and moments of inertia. 

Eq. (^) shows that the deformation of the ring under the action of the force given in Eq. 
(p^), does not depend on twist rigidity a^. Therefore, such external forces do not produce 
any twist and can only lead to bending of the ring. This result remains valid for more general 
distributed forces on the centerline, provided they act only in the plane of the undeformed 
ring. Inserting the expressions for the deviations of the Euler angles, Eqs. ( pBf ) and (0), into 
Eq. (15), we find that 6003 = and consequently the variation of angle ip can be expressed 
in terms of the conformation of the centerline (angle 6), as rd [6iIj{s)] /ds = 66{s). Since the 
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sum of twist and writhe numbers is a topologically conserved number, writhe is invariant 
under such deformations. 

Figure 6 shows the effect of spontaneous (constant) angle of twist, ipo, on the response of a 
ribbonhke (02 ^ ai) ring to compressional forces apphed at opposite points of the centerhne. 
The forces, shown by arrows, he in the plane of the undeformed ring. In the case of a ribbon 
with short axis lying in the plane of the undeformed ring, ipo = 7r/2, the ring remains planar 
in the course of deformation (Fig. 6a). A ribbon with short axis lying normal to the plane 
of the ring, ipo 0, undergoes three dimensional deformation (Fig. 6b). At first sight. 
Fig. 6b appears to suggest that the ring is twisted, in contradiction with the previously 
made statement that its configuration is twist-free. However, as is evident from the Eq. 
(15) for the density of twist, 600^ = d6ip/ds — 66 /r, and from the discussion in Section 
the mathematical definition of the twist of a filament with nonvanishing spontaneous 
curvature (r 7^ 00) involves both the rotation of the cross section about the centerline and 
the curvature of the centerline itself. The fact that the two effects cancel exactly in Figs. 6a 
and 6b is a consequence of the fact that the forces act entirely in the xy plane and do not 
produce a component of torque along the contour of the ring, that could give rise to twist. 

8. Dynamics 

Consider small instantaneous deviations 5x(s, t) = x(s,t) — xo(s) of the centerline of the 
ring from its stress free position, xo(s). We express 5x(s,t) in terms of its projections on 
the triad vectors of the undeformed ring, tofc(s), 

(5x(s,t)=^^5xfc(s,t)tofe(s). (85) 

We proceed to write the Langevin equations that govern the dynamics of fluctuations of 
the centerline, 5x, and the dynamics of angular fluctuations of the cross section about the 
centerline. Sip. Some care should be exercised in deriving the Langevin force from the expres- 
sion for the elastic energy, Eq. (0), since up to this point we have used the inextensibility 
constraint, Eq. (§). In order to avoid complications associated with the introduction of rigid 
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constraints ||2^, we replace the strict inextensibility condition, dSx^/ds = 0002^X1 — uJoiSx2 



(see Appendix ^ , by an energy penalty 



2 



Uext = -^<^ext ds u;o25a;i + uoi6x2j • (86) 

The persistence length aext describes the rigidity of the filament with respect to local com- 
pression and extension. The total elastic energy Utot = U +Uext is the sum of contributions of 
bending and twist modes, Eq. (17), and extensional modes, Eq. (|86[) . We will use the above 
expression for the total energy of an extensible filament in the Langevin equations, and take 
the limit of an inextensible filament {aext/^ ~^ 00) only in the end of the calculation. 
The Langevin equations are 

d'^5ip{s,t) d6ip{s,t) 6Utot , . 

+ '^^T- + 6ms,t)] = ^^^'^ ^^^^ 

Here m and I are mass and moment of inertia (with respect to the centerline) per unit length 
and ? and are translational and angular friction coefficients. The fiuctuation-dissipation 
theorem relates the amplitudes of the random forces f and to these friction coefficients, 

(/,(s,t)) = 0, {f,{s,t)f,is',t')) = 2kBnS,As - s')S{t -t'), (89) 
{^4s,t)) = 0, {^4s,t)^4s',t')) =2kBT<;^6{s - s')6{t -t'). (90) 

In writing the above equations we neglected hydrodynamic interactions and therefore the 
treatment is analogous to the Rouse model of polymer solution dynamics pO |. 



Using the relation between the deviations of the coordinates, (5x, and those of Euler 
angles, 69, 6ip, and 6ip (see Appendix 0), and neglecting rigid body translation and rotation 
of the ring, we rewrite the above Langevin equations in terms of the Fourier components of 
the deviations of Euler angles from their equilibrium values (see Eq. (^l])). 



ar,v{n, t) = + |^(n, t), (91) 



e^(n,t)) = 0, (e^(n,t)ev(-n,t')) = 2kBT6^r,'<ir,Lr,{n)6{t - 1'). (92) 
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The time derivative operators a,, associated with the three Euler angles are, 

ae = = a = m— + and ^ + 

where the corresponding friction coefficients are = = ^ and <^^. The elastic energy that 
appears in the Langevin equations (91), is given in terms of the amplitudes of the Fourier 
modes in Eq. (26). Conveniently, the matrix of kinetic coefficients L is diagonal in the Euler 
angle representation, 

Le{n) = n'^/r^, L^{n) = (n^ - lY/{n'^ + 1^ and L^{n) = (94) 

In the following we proceed to solve the Langevin equations and obtain explicit expres- 
sions for the dynamic correlation function of writhe fluctuations. We focus on this correlator 
since it is an integral characteristic of the ring and is therefore simpler than the two-point 
correlation functions of Euler angles, that depend on the separation between the points. 
Although the general solution of the linear equations can be obtained, we will assume (as it 



is often done in the literature |[T^,^,^) that the relaxation of the twist angle ijj is much 
faster than that of the angles 6 and if. Consequently, we can minimize the energy with 
respect to ij{n,t) and express it in terms of 9{n,t) and (p{n,t). With this substitution, the 
(3 X 3) matrix Q(n) in 6{n, t), <^(n, t), i'{n, t) space, Eq. (^71) , is reduced to a (2 x 2) matrix 
Q'(?T,) in 9{n,t), (p{n,t) space. 



Q'(n) ^ 



( 03^1(^2-1)2 ^303^2(722-1) ^ 



a3n2 + Ai 



(95) 



74303722(722 - 1) (a3y42n2 + 0102)^2 ^ 

As shown in Appendix P, the solutions of the Langevin equations can be expressed in 
terms of the eigenvalues Ai 2(^) of the matrix P(n) of the linear form 

^vi^'^j-ip;;^ = ^ri'^vv'i^Win,t), r],r]' = 9,^ (96) 

that appears in the Langevin equation, (91). Explicit expressions for these eigenvalues are 
given in Appendix p, Eq. (B2). 
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Taking the Fourier transform of Eq. (p9| ) (Appendix P) with respect to the frequency 
a;, we find the two-time correlation function of the Fourier components of Euler angles, 

1 



{ri{n,t)ri (-n,0)) 



^^^^^ _ {[Mn)9iin,t) - Ai{n)g2in,t)] x 

{f]{n)fi'{-n)) -kBTLn{n) [gi{n,t) - g2in,t)]6nn'} , 



(97) 



where gk{n,t) describes temporal decay of correlations of normal modes with wavevector 
27m/r {gk{n,0) = 1), and where {fj{n)fj' {—n)) is the previously calculated equal-time equi- 
librium correlation function (see section 

In the inertial limit (i.e., for modes with inertial time scale shorter than viscous relaxation 
time), 4mAfe(n) > <j^, the function gkiji^t) describes damped oscillations with characteristic 
frequency uj^iji), 



cos ujk{n)t + - — — 
2muJk{n) 



sin Uk{n)t 



exp ( — 



2m 



m 



Am? 



(98) 
(99) 



The characteristic relaxation time is 2m/ independent of the wavelength of the mode. At 
short times, t <^ m/<^. 



gk{n,t) ~ cos 



\jKk{n)/mt 



(100) 



In the dissipative limit, 4mAfc(n) < the function gk{n,t) describes pure decay of 
correlations, with characteristic times Ti(n) and T2(n), 

t 



9k{n,t) 



n I t 

exp 



n - T2 



T2 



exp 

n / T-2 - Ti \ T2 



± 



m 



(101) 
(102) 



2A,(n) \AAlin) Afc(n)- 

In the limit of negligible inertia, m 0, we get T2{n) — > and the relaxation can be 
described by simple exponential decay. 



gk{n,t) ~ exp [-tAk{n)/q] 



(103) 
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The dynamic writhe correlation function is derived in Appendix ^ 

{5Wr{t)5Wr{0)) = Wr\n,t)= ^ Wr\n)gi{n,t)g2{n,t) (104) 

n^O,±l n^O,±l 

where the equihbrium mean squared amphtude of fluctuations of writhe, Wr'^{n), is given in 
Eq. (|6^ . In the limit of negligible inertia, m — > 0, we substitute Eq. (|103|) for gk{n,t) into 



Eq. (|104|) , and find that the writhe correlation function for mode n decays exponentially 



with time 

Wr'^{n, t) = Wr^{n)e-^^^^''\ (105) 

The characteristic decay time T{n) is given by 

r \ n^ + l a^n^ + Ai 
T[n) = — — — -r- . (lOo) 

ttkbT n^[n^ — 1)^ (ai + a2)a3n'' + Aia^ + aia2 

In the short wavelength limit, n ^ 1, only the sum of bending persistence lengths, Oi + 02, 
appears in T{n). Indeed, on small scales the filament behaves as a straight inextensible rod 
whose properties do not depend on the twist persistence length 03, or on the spontaneous 
twist angle ifjQ. 

In Figs. 7 and 8 we plot the writhe correlation function as a function of time measured 
in units of (jr^/vr/csT, in the inertial regime, for 2'KmkBT /q^r"^ = 10. Its Fourier transform 
plotted as a function of the frequency u measured in units of TcksT/qr'^ , is shown in the 
insert, on the upper right side of the figure. In Fig. 7 the parameters correspond to a 
circular cross section and identical persistence lengths, di — (I2 — 

= 2r. Oscillatory decay 

of writhe correlations as a function of time is observed, but the correlation remains always 
positive. A small number of fundamental frequencies can be detected in the oscillatory 
pattern, and identified with peaks observed in the frequency spectrum. The amplitudes of 
these peaks decrease monotonically with the frequency, and the largest peak is at = 0. 
The case of asymmetric cross section and dominant bending rigidity, ai = = 2r and 
02 = 20r {ipQ = 7r/4) is shown in Fig. 8. The correlation function decays rapidly to zero 
and, at later times, oscillates between positive and negative values. Since for t <^ m/<^, 
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dissipation is negligible, the fast initial decay of correlations is the result of dephasing of the 
oscillatory contributions of a large number of modes. In the frequency domain, there is no 
peak at = and the peak amplitudes have nonmonotonic dependence on the frequency. 

In Fig. 9 we plot the writhe correlation function in the dissipative regime where inertia 
is negligible, 2'KmkBT /q^r"^ = 10~^, for ribbonlike rings with different bending and twist 
rigidities. The amplitude of writhe fluctuations is smallest for a ribbon whose shorter axis 
of inertia is normal to the plane of the ring (see Fig. 6b). The amplitude increases by more 
than a factor of two when the shorter axis of inertia lies in the plane of the ring (see Fig 
6a). Twist rigidity decreases the fluctuation amplitude but the effect is rather weak. Since 
inertial oscillations are completely suppressed in this overdamped regime, the correlations 
decay monotonically with time. The curves exhibit fast short time relaxation, followed 
by an exponential decay, in qualitative agreement with numerical simulations reported in 
reference ||19|]. An analytic expression for the time dependence of the correlator at short times 
can be derived (see Eq. ( [BT8D in Appendix [§). The predicted (^[6Wr{t) - 5Wr{0)f) oc t^/^ 
dependence of the writhe correlation function and the fact that it depends only on the 
bending rigidity (ai and 02), are direct consequences of the observation that at short times, 
the relaxation is dominated by straight rod contributions to the spectrum (r^q) oc l/q'^)- 

The above results can be directly applied to the study of deformation of macroscopic 
rings by external forces and torques. Unlike the case of microscopic filaments where dis- 
sipation dominates inertia and only overdamped behavior is expected, inertial effects play 
an important role in the dynamic response of macroscopic objects (for this reason they 
were included in the preceding analysis). According to the fluctuation-dissipation theorem, 
dynamic correlation functions of Euler angles can be treated as generalized susceptibili- 
ties that determine the response of the ring to externally applied torques and forces. The 
time-dependent generalization of the static relation between deformation (in terms of the 
deviation of the Euler angles from their equihbrium values) and applied force, Eq. (^), is 

Sr^{s,t) = ^ £jt'Y.^, f ds'^ {Sv{s,t)6r^\s',t')) M,,{s',t'), (107) 
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where the moments M due to external forces, are given by Eq. (80). The relaxation of the 
deformation following the release of external moments at time t — 0, M(s, t) — M(s)6'(— t), 
for t > is given by 



9. Discussion 

In this paper we presented the statistical mechanics of fluctuating rings. We derived analyt- 
ical expressions for various static properties of such rings, including two-point correlation 
functions of Euler angles, correlation function of tangents to the ring, rms distance be- 
tween points on the ring contour, radius of gyration and probability distribution function of 
writhe, as function of persistence lengths associated with bending and twist deformations of 
the ring. We found that the amplitudes of fluctuations of the Euler angles and 9 diverge 
in the limit of vanishing twist rigidity. We would like to emphasize that the situation differs 
from the case of straight filaments for which the twist density duj^^'^ = dSip / ds depends only 
on the fluctuations of the Euler angle ■0- For such filaments, vanishing twist rigidity (03 = 0) 
implies that there is no energy penalty for twisting the cross section about the centerline, 
but the presence of bending rigidity (ai,a2 7^ 0) suffices to suppress spatial fluctuations of 
the center fine about its straight stress free configuration. Thus, if we are only interested 
in the statistical mechanics of the spatial conformations of the centerline, accounting for 
bending rigidity suffices to provide an accurate description of straight fluctuating fllaments. 
For rings, inspection of the elastic energy, Eq. (17), shows that when 03 = fluctuations 
with d59/ds + 5ip/r — have zero energy cost and, since in the absence of twist rigidity 
the angle 6ip can always adjust itself to satisfy the condition 6ip = —r dS9/ds, there is no 
elastic energy penalty for out-of-plane fluctuations of the ring and the amplitude of such 
fluctuations diverges even if the bending rigidity remains flnite. Therefore, wormlike chain 
theories in which only bending rigidity is taken into account, can not model the spatial 
conformation of fluctuating rings. 




(108) 
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We found that a crossover length scale — f\Jaz/ah exists, below which straight rod 
behavior dominates and writhe of the centerline and twist of the cross section about it are 
decoupled, and above which spontaneous curvature becomes important and twist affects 
the three dimensional configurations of the centerline of the ring. In this context we would 
like to propose the as yet unproven but plausible conjecture, that the existence of this 
crossover does not depend on the topology of the ring and is characteristic of filaments with 
spontaneous curvature in their stress free state. 

Although the main focus of this work is on the statistical mechanics of fluctuating rings, 
we used the fluctuation-dissipation theorem in order to predict mechanical response to ex- 
ternal torques and forces, and showed that the deformation of ribbonlike rings depends in 
an essential way on the orientation of the cross section in the stress free state. Finally, 
we derived the Langevin equations that govern the dynamics of fluctuating rings, and cal- 
culated the two-time correlation function of writhe fluctuations. Depending on the values 
of the parameters, one can move from an inertial regime where relaxation is accompanied 
by temporal oscillations, to a non-oscillatory, purely dissipative regime. In the dissipation 
dominated range, the relaxation at short times is determined by bending rigidity only. This 
agrees with the expectation that short time relaxation is dominated by small scale, straight 
rod behavior. At longer times, the decay of writhe correlations depends on both bending 
and twist rigidities. While inertial effects are not expected to be important for microscopic 
objects such as small plasmids, our dynamic response functions can describe the relaxation 
of rings of arbitrary mass and size, following the cessation of externally applied forces and 
torques. 

We would like to comment on the limitations of the approach presented in this paper. 
The domain of applicability of our theory is limited to the weak fluctuation regime, in the 
sense that the deviations of the Euler angles from their values in the undeformed ring, must 
be sufficiently small. Although, in principle, our general formalism is applicable to rings with 
arbitrary spontaneous twist in their stress free reference state, the analysis of this problem 
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meets with considerable mathematical difficulties and is the subject of ongoing work. Finally, 
we would like to emphasize that since our theory is based on the linear theory of elasticity of 
thin rods, all persistence lengths and radii of curvature are assumed to be much larger than 
the diameter of the filament that serves as a small scale cutoff. Consideration of microscopic 
physics on length scales smaller than this diameter requires the introduction of additional 
model assumptions (see references [p3H25[|) and is beyond the scope of this paper. 

After this work was submitted for publication we learned about a related study of thermal 
fluctuations in DNA plasmids in which the writhe distribution function was also calculated 
||26| . This work is complementary to ours: while we assume that the equilibrium stress free 
state of our filament is that of a planar untwisted circular ring, reference deals with 
filaments with straight untwisted stress free state. Conceptually, the ring is then formed by 
bringing the ends together and sealing them, with or without the addition of twist. Since 
such rings have locked-in internal stresses, the two procedures are nonequivalent in general. 
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Appendix A: Derivation of Langevin equations 

For small deviations from the stress free state, the variation of triad of vectors can be written 

as 

6ti = - {66 cos ipo + 5ip sin ipo) tos + 6ipto2, 

5t2 = {69 sinipo - 5(p cos ipo) tos - 6i>to2-, (Al) 
(5t3 = — {56smtpQ — 5ip cos tpo) to2 + {56cosipQ + 5ipsm.ipQ) toi, 

where the vectors toi are defined in Eq. (|T^). Substituting Eqs. ( [AID i^^o ^^e inextensibility 
condition, Eq. (|^), we obtain the following equations for the deviation 5x(s) of the position 
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vector of a point s on the ring contour, from its value in the undeformed state 
d6x3 



ds 
d6x2 



- UJo25xi + Co'oi(5x2 = 0, (A2) 

+ UJ016X3 — UJQ36X1 = 66 sin ipQ — 6ip COS ipo, (A3) 
as 

+ 10026x3 - UJ036X2 = 66 cos i/jQ + 6ip sin iI^q. (A4) 



Eq. ( [A^ ) is the hnearized form of the inextensibihty condition, and Eqs. ( [A3D and ( |A4| ) 
relate the deviations of Euler angles to those of spatial positions. 

Fourier transforming Eqs. ( |57D and (^^), yields the Langevin equations for the Fourier 
components x (n) and '?/'(n), 

6U 

Q;x(n, t) + —, — + /i(n, t)c*{n) = f (n, t), (A5) 

a^^{n,t) + -j- - = ^^{n,t), (A6) 

dip[—n, t) 



where a and are defined in Eq. (|93|) and 

t) = ^ — x(r;,, t) ■ c(n), c(n) = ( — sin?/'o, — cosV'o, ^"-)- (A7) 

In the limit aext — cxd, ^{n, t) can be considered as a Lagrange multiplier that accounts for 
the inextensibihty condition, Eq. ( |A2| ) or, equivalently, for its Fourier transform, x(n, t) ■ 
c(n) = 0. The correlators of random forces in Eqs. (A5) and (A6), take the form, 

= 0, = 2^T5,,5(t-t'), (A8) 

{i^{n,t)) = 0, (f^(n,t)e^(-n,t')) = 2c;^T5(t - 1'). (A9) 



Using Fourier transforms of Eqs. ( |A^ ) - (^J) we rewrite the Langevin equations (A5) - 
(A9) in terms of Euler angles, Eqs. (91) and (92). 
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Appendix B: Solution of Langevin equations 



In order to find the solution of the Langevin equations, we first calculate the eigenvalues 
and eigenfunctions of the matrix (see Eq. (PO) for its definition) , 



^ ^303 ^"^2^^ {,A2a^n^ + 0102) ^^Tipi y 



(Bl) 



Eigenvalues of this matrix have the form 



-KksT n^{n^ - 1)^ (ai + a2)a^n^ + Aia^, + 0102 ± A 
^^'^^"^ = 2r3 + 1 ^^;I^T:4^ ' ^^2) 

= [(fli — 02)03^^ + ^103 — 0102]^ + 4n^a2a3(ai — a2)(oi — 03) sin^ 

The eigenvalues A/(.(n) vanish when n = 0, 1. These modes are associated with rigid body 
rotations of the ring and are not considered further below. In the limit \n\ ^ 1 both eigen- 
values increase with the fourth power of n, Ak{n) ~ irkBTakn^/r^. This dependence of 
the eigenvalues (g = 2im/r is the wavevector corresponding to the n— th mode) is character- 
istic of bending fiuctuations of straight rods, in accord with the expectation that small-scale 
fiuctuations of a ring are indistinguishable from those of a straight rod. 

The matrix P(n) becomes diagonal (and the angles 9 and if become decoupled), in the 
case of a circularly symmetric cross section (^3 = 0), when t/'o = or 7r/2, and in the limit 
03 — > 0. In all of these cases the mode Ai(n) describes both fiuctuations perpendicular to 
the plane of the ring and twist fiuctuations, and the mode A2(n) describes fiuctuations in 
the plane of the ring. Since Ai(r2) when —>■ 0, twist fiuctuations destroy the circular 
shape of the ring in the wormlike chain model, where twist rigidity is not taken into account. 

The eigenvalues take a particularly simple form for ipo = 0, 

TrksTn'^in^ -Ifaitts . , . TikBTn^{n^-lfa2 
Ai n = — I , A2 n = — § • B3 



Since the matrix P(ti), Eq. (|B1|), is not symmetric, it has different right and left 
eigenfunctions. We denote the right eigenfunctions, corresponding to eigenvalues Afc(n), 
by fjkin) = {9k{n),(pk{n)}, (where k = 1,2). Left eigenfunctions can be written as 
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1k{~n)L^^{n). Since the matrix P(n) is real, we have rik{—n) = rjKn). For each n 7^ 0, ±1, 
the above eigenfunctions are normahzed by conditions: 

J2^L-\n)r]k{n)r]i{-n) =6ki, ^(n)/]^ (-n) = L^(n)5^^/. (B4) 

Expanding the matrix P(?7,) over its eigenfunctions we find: 

The solution of Eq. (91) can be found by Fourier transforming it with respect to the 
time t, and substituting Eqs. (p6|) and (|B5|). This yields 



CtuJ = —mcu'^ + KLO, (B7) 

where the correlators of the Gaussian random force are 



Calculating the correlation functions of Euler angles, we get 

where rj, rj' = 6, ip. In the time domain this gives the following expression for the dynamic 
correlation functions of the Fourier transforms of Euler angles, 

{f,{n,t)v'{-n,0)) = fci.r5:^^^^^^|pl^,(n,t), (BIO) 

where 

2qAk{n) du cos ujt 
9k{n,t) = / — — (Bll) 

The function gk{n,t) describes temporal decay of correlations of normal modes, with 
wavevector 27rn/r. One can verify that the for t = 0, integration gives gk{n,0) = 1, in- 
dependent of m and Instead of calculating the eigenfunctions fik{n), we notice that the 
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combinations f]k{n)f]'i.{—n) / Ak{n) in the above expression, can be evaluated using the pre- 
viously derived expressions for the equilibrium (equal time) correlators {fj{n)fj' {—n)) , and 
the normalization conditions, Eq. (|B4|) . Substituting these expressions into Eq. ( |B1CI| ), we 
arrive at Eq. (97). 

We now turn to writhe fluctuations (see Eq. (|5P|)), 

SWr{t) = -T.J^ J ^^~-'^-^)^^^'*' ^^12) 

and proceed to calculate the dynamic correlation function of these fluctuations, 

{6Wr{t)6Wr{0)) = J2 Wr\n,t), (B13) 

where the contribution of mode n is 

Wr\n,t) = k%T^n^ [ ^ cos[iuj + uj')t]x 

J Zn J ZTT 

{0^{-n)9^^{n)^ {(fcu'{n)(f-uj'i-n)) - (^9^{-n)(p-^{n)^ {o^uii{n)(p^i{-n)^ . ^^^^^ 
Substituting Eq. (43) for the correlation functions of Euler angles yields, 

Ak{n}Ak'{n} 

9k{n)9k{-n)ipk'{n)^k'{-n) - ek{n)^k{~n)ipk'{n)Ok'{-n) . 

(B15) 

The only nonvanishing contributions to the sum in Eq. (B15) are k = l,k' = 2 and 
k = 2, k' = 1, and both have the same value. As a result, we flnd that the contribution of 
the n— th mode to the dynamic correlation function can be recast in the form 

Wr^{n, t) = Wr^{n)gi{n, t)g2{n, t), (B16) 

where Wr'^{n) is the mean squared amplitude of the n— th mode of writhe fluctuations in 
equilibrium, calculated earlier in Eq. (|62D. 

Finally, we would like to comment on the short time behavior of the correlation function 
in the dissipative regime. Using Eq. ( p.05| ) we flnd 



mm - mm?) - ^ I ^ [1 - e—'] . (BIT, 
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At small times, t <^ t(2), this sum is dominated by terms with n ^ 1. Replacing the sum 
by an integral we find 



[SWr{t) - SWr{0)f 



4^2 



3-' 



ai + a2)t 



1/4 



(B18) 



7r^aia2 

where F is the gamma function. The characteristic relaxation rate 7rkBT{ai+a2) / <ir^ depends 
only on the bending rigidity of the ring. 
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Figure captions 

Fig. 1: Plots of two-point correlation functions of Euler angles {6ri{s)6ri'{0)) vs. the 
contour distance between the points, s, in the interval < s < 27rr: {59{s)S6{0)) (cross), 
{Sip{s)Sip{0)) (diamond), (5V'(s)5V'(0)) (circle) and {59{s)Sip{0)) (solid line). The parameters 
arc ipo = and ai = 02 = lOr, as = lOOr. 

Fig. 2: Plots of two-point correlation functions of Euler angles {5r){s)Sr)'{0)) vs. the 
contour distance between the points, s, in the interval < s < 27rr: {S9{s)S9{0)) (cross), 
{dip{s)6ip{0)) {circle) and {69{s)5ip{0)) (solid line). The parameters are -^0 = and ai = 
a2 = lOr, as = r. 

Fig. 3: Plots of nondiagonal two-point correlation functions of Euler angles 

{6ri{s)6ri'{0)) vs. the contour distance between the points, s, in the interval < s < 27rr: 
{59{s)5(p{0)) (cross), {5(p{s)5i/j{0)) (circle) and (5^ (s) (5-0(0)) (solid line). The parameters are 
■00 = 7''/4 and oi = lOr, 02 = lOOr, 03 = lOr. 

Fig. 4: Plot of dimensionless rms distance between points on the ring contour 
[x(s) — x(0)]^^ /r^ vs. the contour distance between the points, s, in the interval < 
s < 27rr. The parameters are ipo — and: ai — a2 — — lOr (solid line), ai — a2 — lOr, 
as = r (box), ai = a2 = r, as = lOr (cross) and ai = lOr, a2 = as = r (diamond). 

Fig. 5: Plot of probability distribution function of writhe p{x, Ai/as) vs. x — 
(TTy^oIo^/r) 5Wr for ^1/03 = 0.1 (solid line), 5 (cross) and 20 (box). 

Fig. 6: Plots of deformation of ribbonlike rings (with a2/ai = 10"^) by compressional 
forces (see arrows), a) '4>q — t^/'2. and b) -^o — 10~^. 

Fig. 7: Plot of dynamic correlation function of writhe fluctuations {5Wr{t)5Wr{Q)) vs. 
time t (in units of qr"^ /TrkBT), for a ring with a circular cross section and persistence lengths 
ai — a2 — as — 2r, in the inertial range 27r'mkBT / q^r^ — 10. Plot of the Fourier transform 
of the correlation functions vs. frequency cu (in units of tt/cbT/^t^) is shown as an insert in 
the upper right hand side of the figure. 

Fig. 8. Plot of dynamic correlation function of writhe fluctuations {SWr{t)dWr{0)) vs. 
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time t (in units of qr^ /irkBT), for a ribbonlike ring with persistence lengths oi = 03 = 2r 
and 02 = 20r {i/jq = 7r/4), in the inertial range 27rmA;sT'/^^r^ = 10. Plot of the Fourier 
transform of the correlation functions vs. frequency lu (in units of TrkBT/qr^) is shown as 
an insert in the upper right hand side of the figure. 

Fig. 9. Plot of dynamic correlations function of writhe fluctuations {5Wr{t)5Wr{0)) vs. 
time t (in units of qr^/frkBT), for ribbonlike rings in the dissipative range — 
10~^: The parameters are t/jo — and oi = 2r, 02 = 20r, 03 = 2r (cross), oi = 20r, 
02 = 2r, 03 = 2r (circle) and oi = 20r, 02 = 2r, 03 = lOr (diamond). 
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